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For any near-threshold asymptotic regime and for any Feynman diagram (involving loop and/or 
phase space integrals), a systematic prescription for explicitly constructing all-logs, all-powers 
(all-twists) expansions in perfectly factorized form with explicit integrals for coefficients, is pre- 
sented. The distribution-theoretic nature of the method of asymptotic operation employed allows 
treatment of totally exclusive phase space in the same manner as loop integrals. 



Introduction 1 

The purpose of this Letter is to present a summary of sys- 
tematic recipe for construction of asymptotic expansions of 
Feynman diagrams near threshold values of kinematic vari- 
ables. The recipe lays a foundation, heretofore lacking, for sys- 
tematic higher-order calculations as well as for all-order con- 
struction of asymptotic expansions in operator/effective La- 
grangian form. The range of applicability of the recipe includes 
many concrete applications such as near-threshold production 
of, say, electron-positron pair in QED (and many similar prob- 
lems in the Standard Model); the small-x problem in deeply 
inelastic scattering; etc. 

The recipe is a culmination of the development of the tech- 
nique of asymptotic operation (A.y-operation or simply AO) 
which has been the driving force behind the continuous prog- 
ress in the systematic 1 studies of asymptotic expansions of 
Feynman diagrams in masses and momenta since before 1982 
[l] 2 . The Euclidean variant of AO yielded powerful calcula- 
tional formulas for the short-distance expansion in the MS 
scheme [3; 4] and for mass expansions [5-8] that were used in 
a number of NNLO calculations in QCD and the Standard 
Model (e.g. [9-11]). The non-Euclidean extension of AO pre- 
sented below is intended to play the same role with respect to 
the near-threshold problems as the Euclidean variant did with 
respect to those calculations. 

Roughly speaking, Euclidean regimes correspond to cases 
when some masses and momenta are componentwise larger 
than others. 3 On the other hand, one deals with a truly non- 
Euclidean regime when some momenta have both large and 
small components, and the large components approach special 
values that result in a non-trivial overlap of singularities of dif- 
ferent factors in the momentum-space integrand. Such situa- 
tions correspond to thresholds from the point of view of analy- 



1 "Systematic" here means all-logs, all-powers treatments complete with 
explicit calculational formulas in a form maximally suitable for practical 
applications. 

2 For a review and complete references to the original publications, as well 
as for comparison with the conventional methods such as the BPHZ theory 
and the technique of leading logarithmic approximation, see [2]. 

3 To this class also belong degenerate cases with external large momenta 
are fixed on mass shell which simply corresponds to complex-valued large 
external momenta from the point of view of Euclidean space. 



ticity properties of Feynman diagrams; cf. the discussion 
in [12]. 

The central difficulty of constructing asymptotic expansions 
of Feynman diagrams in masses and momenta is that formal 
Taylor expansions of integrands possess non-integrable singu- 
larities localized on variously intersecting manifolds in the 
space of integration momenta. The key observation (from 
which the entire theory of AO unfolds in a logical manner) is 
that the difficulty is a manifestation of the distribution- 
theoretic nature of the expansion problem [l] 4 , and that the 
crucial mathematical task is to find expansion of the integrand 
in the sense of distributions 5 . Asymptotic operation is a pre- 
scription that yields such an expansion for a given integrand (= 
product of propagators). 

The structure of AO is fixed by the extension principle 
[1; 15] — a very general but essentially simple proposition re- 
lated to the well-known Hahn-Banach theorem (see e.g. [17]). 
The resulting prescription is, roughly, as follows (see also 
Sec. 2 below): 

(i) The formal expansion of the integrand should be supple- 
mented with counterterms — linear combinations of 
(derivatives of) <5-functions with coefficients that depend on 
the expansion parameter (cf. below Eq.2.9). 

(ii) Concrete integral expressions for the coefficients valid 
within the precision of expansion are obtained from the so- 
called consistency conditions [18; 15]. 

(iii) A fine-tuning of the consistency conditions to achieve 
purely power-and-log dependence of the coefficients of coun- 
terterms on the expansion parameter (the property of perfect 
factorization [3]). 6 

After the counterterms of AO are found, obtaining the corre- 
sponding all-order operator expansions is a matter of more or 
less straightforward combinatorics (cf. the Euclidean case [8]). 



4 This observation was influenced by Bogoliubov's analysis of the UV 
problem [13; 14]. 

3 For precise definitions see [15; 16]. 

6 Note that the deterministic logic of AO leads one step-by-step towards the 
solution. Contrast this with the BPHZ-type methods where one, in fact, 
needs to guess the result (the forest formula) as a whole in order to proceed 
to, say, formal proofs. For instance, the BPHZ-type interpretation of the 
general Euclidean expansions needed to borrow the explicit results that 
had been previously obtained within the framework of AO (for details 
see [19]). 
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Thus, the analytical focus of the entire theory of asymptotic 
expansions of Feynman diagrams is in finding the coefficients 
of the counterterms of AO in a form best suited for applica- 
tions. The point is that, in general, asymptotic expansions are 
not unique (which is reflected in an arbitrariness of the expres- 
sions for coefficients of AO obtained from consistency condi- 
tions). Uniqueness is restored [15], however, if one imposes on 
expansions the requirement of perfect factorization [3] — 
which at the level of individual diagrams stipulates that the 
expansions run in pure powers and logarithms of the expansion 
parameter. Apart of its importance for applications 
(determination of power-suppressed corrections is impossible 
without it [20]), the uniqueness of power-and-log expansions 
has a number of technical advantages: 

(i) AO commutes with multiplications by polynomials thus 
allowing one to ignore complications due to non-scalar 
particles; 

(ii) expansions inherit all algebraic properties of the initial 
integrands (such as the Ward-Takahashi-Slavnov-Taylor iden- 
tities); 

(iii) maximal calculational simplifications; (iv) considerable 
advantages for a regularization-independent treatment [16; 21], 
which is potentially important for supersymmetric models. 

The consistency conditions have the form of integrals in- 
volving some of the propagators of the original Feynman inte- 
grand as well as a test function whose behavior at the singular 
point is fixed but which is otherwise arbitrary (cf. Eq. 2.12). 
For the Euclidean case, it was found [18; 15] that it is suffi- 
cient to replace such test functions with suitable polynomials 
(to be understood as an appropriate limiting procedure; for a 
rigorous treatment see [16; 21]). It turned out that such a re- 
placement modifies the expressions only within the precision 
of the expansion (which is always allowed), whereas the sim- 
ple scaling of the integrand in momenta and masses ensures 
the power-and-log dependence of the resulting expression on 
the expansion parameter. This allowed us to include into sys- 
tematic consideration the entire class of Euclidean asymptotic 
regimes [18] and subsequently to expand the scope of operator- 
expansion methods to such regimes. 

As was emphasized from the very beginning [18; 15; 2], the 
scenario of AO is completely general and by no means limited 
to Euclidean cases. The specifics of the non-Euclidean regimes 
is that the singularities of integrands are localized on non- 
linear manifolds (light cones and mass shells) and that non- 
zero finite limiting values for external momenta break the 
usual scale invariance of integrands. As a consequence, the 
simple trick that yielded power-and-log dependence in the 
Euclidean case is no longer sufficient. However, it is important 
to understand that the consistency conditions are obtained from 
first principles without any restricting assumptions, and there- 
fore possess all the flexibility to accommodate any additional 
requirement that one may lawfully impose — in particular, to 
perform an appropriate fine-tuning to achieve the required 
power-and-log dependence on the expansion parameter. 

The resulting problem and its solution were identified in 
[2]: the problem consists in occurrence of the so-called oscu- 
lating singularities' ', whereas the solution is given by the so- 



7 i.e. singularities whose singular manifolds touch rather than intersect in a 
general fashion; for a detailed discussion see [12]. When a mixture of os- 
culating and transverse intersections occurs, simple uniform scaling rules 
for description of the strength of the singularity no longer suffice. 



called homogenization — a secondary expansion that comple- 
ments the consistency conditions for the coefficients of AO. 
The homogenization splits those coefficients into pieces with 
strict power-and-log dependence on the expansion parameter 
(powers needing not be integer). Below we expand the sce- 
nario of [2] by presenting explicit rules for the homogeniza- 
tion. 8 

It is remarkable that a self-contained analytical recipe for an 
individual diagram can be summarized in a rather compact 
universal form for arbitrary non-Euclidean regimes despite 
their larger analytical variety than in the Euclidean case. Such 
a description would not be possible without using the language 
of AO. However, there is both an increasing familiarity among 
physicists with the technique of AO, and a tendency to use it in 
various non-standard physical problems due to its power and 
flexibility (cf. [22-26]) 9 . On the other hand, not all such works 
have been equally successful. 10 Therefore, it would be useful 
— prior to a more complete treatment which would require a 
substantially longer text — to give a summary of the procedure 
in a form suitable for calculations that can also serve as a 
starting point for derivation of all-order operator-form expan- 
sions for various regimes. 

Before we proceed to formulas, a few remarks are in order. 

(i) One sometimes uses the term Minkowski space regimes 
(cf. [27]). However, Minkowski space per se allows both 
Euclidean regimes (cf. their treatment in Minkowski space in 
[12]), and a more complex class of non-Euclidean, or near- 
threshold regimes. I emphasize the distinction in order to avoid 
confusion due to vagueness of terminology and argumentation 
in many publications on the subject. 

(ii) In some recent publications the term "threshold expan- 
sion" was misused to denote cases with kinematic parameters 
set exactly at threshold values with expansions running with 
respect to some other parameter (e.g. internal mass), or degen- 
erate thresholds tractable by Euclidean methods. In this Letter 
we consider true non-degenerate thresholds (including non- 
zero ones) intractable by Euclidean methods, and expansions in 
a parameter that measures closeness to such a threshold. 

(iii) The method of AO considers integrands in momentum 
space as distributions prior to integration. It is therefore ideally 
suited for studying physical problems where integration over 
the phase space of final state particles should not be performed 
in an explicit fashion. Indeed, the 5-functions that describe the 
phase space — e.g. 6(p Q )8(p 2 - m 2 ) etc. — are, from the 
distribution-theoretic point of view, equally acceptable factors 
alongside the standard causal propagators (p 2 - m 2 + i0) _1 . 
This opens a prospect for a systematic treatment of e.g. the 
problem of power corrections in jet-related shape observables 
in the context of precision measurements of a s according to 

the theoretical scenario outlined in [28] in analogy with the 
case of total cross section of e + e~ — > hadrons where the op- 

8 Ref.[2] outlined the worst-case scenario for non-Euclidean AO because a 
universal description of its structure was not yet available. Thus, the com- 
binatorial complexity due to the homogenization seems now to be less se- 
vere — at least for some asymptotic regimes — than anticipated then. 

9 The integral version of AO [6; 8] (cf. also its first special case — the 
formulas for OPE coefficient functions in the MS scheme [3; 4]) is useful 
in situations where efficient automation is necessary, such as higher-order 
calculations etc. 

10 For instance, the description of AO in ref.[25] is incorrect whereas the 
authors of [26] independently found via AO a correct treatment of a con- 
crete near-threshold problem. 
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erator-product expansion can be used for that purpose [20], 
Problems with exclusive phase space occur on a massive scale 
in physical applications." The various cuts used for event se- 
lection are equivalent to various weights in phase space inte- 
grals, which means that the corresponding matrix elements 
squared are effectively treated as distributions in the momenta 
of the final state particles — a perfect setting for application of 
the distribution-theoretic technique of AO. 

Description of the method 2 

The prescriptions of AO are best described as a formal deri- 
vation rather than a final formula or theorem. The reason is 
that every step of the derivation has a simple concrete meaning 
enabling one to control correctness of formulas in each con- 
crete situation, whereas blindly using a cumbersome final for- 
mula may result in gross errors. 12 

We follow the notations of [15; 16] and focus here only on 
the most difficult — analytical — aspect of the expansion 
problem; the diagrammatic interpretation depends on a con- 
crete asymptotic regime and is a much simpler (combinatorial) 
issue anyway. We describe AO in the form with an intermedi- 
ate regularization similarly to how the Euclidean case was 
treated in [15]. It is not difficult to convert the formulas into a 
regularization-independent form similar to [16] (further details 
specific to the non-Euclidean are given in [30]). In general, we 
have in view a combination of dimensional [31] and analytical 
[32] regularizations; the latter may be needed in the cases 
when the dimensional regularization alone is insufficient (cf. 
the example of [23]). 

The expression to be expanded 2.1 

The collection of all integration momenta (loop and phase 
space) is denoted as p . The integrand of the diagram to be ex- 
panded in a small parameter K is represented as follows: 

G ^ = I\ geG ^' g(P) = ^g{lg(P>K)), 2.2 
where A g (z) = (z±i0) _1 or 8(z) (each such factor will be re- 
ferred to as "propagator"); l„(p,K) is a second order polyno- 
mial of the momenta p and masses (first order polynomials are 
allowed as a special case; cf. gauge boson propagators in non- 
covariant gauges). Note that we allow as factors causal propa- 
gators, their complex conjugates, and phase-space 5-functions 
(for simplicity, the ©-functions of phase space factors are 
omitted; to include them, it is sufficient to modify A accord- 
ingly in all formulas; cf. the example in Sec. 3.1). Various 
polynomials that may occur in the numerator (due to non-scalar 
particles and interactions with derivatives) are ignored because 
AO commutes with multiplication by polynomials [15]. The 
expression 2.2 comprises as special cases loop and unitarity 
diagrams. 



11 1 am indebted to S. Jadach for explaining to me this point. 

12 Incidentally, Collins et al. [25] attempted to describe a non-Euclidean 
example of AO. The result is a bizarre text which makes very little sense 
beyond vaguely echoing the discussions I had with John Collins during my 
three visits to Penn State (I believed we were discussing application of the 
technique of AO to the Sudakov problem; cf. [23]). Several of the formu- 
las of [25] (together with the accompanying textual descriptions) are sim- 
ply incorrect. This calls for a critical reexamination of the "proofs" of 
QCD factorization theorems, which I intend to do elsewhere [29]. 



UV renormalization is assumed to be performed in 
a massless scheme of the MS type and is treated following 
[33; 34; 21] as a subtraction from momentum space integrand 
of its asymptotic terms (in the sense of distributions) for 
p-^oo (see [15; 16] for an exact interpretation of the large-/? 
limit involved). For practical purposes, it is sufficient to 
employ the MS scheme [35] (or any of the massless renormali- 
zation schemes), and treat unrenormalized UV-divergent 
integrals formally as convergent in the usual fashion (cf. the 
prescriptions of the Euclidean AO [15; 8]; a rigorous treatment 
of why this is possible is given in [16; 21]). 

Expansion parameter and asymptotic regimes 2.3 

The expression 2.2 depends on external parameters such as 
masses, momenta of incoming particles etc. It is assumed that 
some of the momentum components and/or masses are small 
compared to others. The small parameter (one with respect to 
which the expansion is to be performed) will be denoted as K. 
In general, one assumes that some of the external parameters 
— masses or momenta — tend to specific values (zero or not), 
and that the differences between the external parameters and 
their limiting values are of order K (extension to cases with 
several scales of the form K n is straightforward). The limiting 
values of external momenta need not be zero componentwise. 
Some examples of the constructs mentioned in the description 
given below, are presented in Sec. 3 (for further examples see 
[23] and [36]). 

Our prescriptions are valid irrespective of what kind of 
threshold the chosen asymptotic regime corresponds to — per- 
haps, none at all in which case there will simply be no singu- 
larities requiring addition of non-trivial counterterms. 

Note also the following rule for the problems with explicit 
phase space: If some phase space momentum components are 
to be treated as small, i.e. 0(k) , they should be made O(l) by 
appropriate rescaling before applying the procedures of AO. 

Formal expansion 2.4 

The construction of AO begins with the formal (usually, but 
not necessarily, Taylor) expansion of the integrand in powers 
of k: G(p) — > 7 K ° G(p) . Each factor is expanded separately, 
the results are formally multiplied and reordered in increasing 
powers of K. The terms of the resulting series possess, in gen- 
eral, non-integrable singularities which have to be examined in 
the geometrical and analytical aspects. 

Geometric classification of singularities 

of the formal expansion 2.5 

Each formal expansion T K °g(p) of each factor g(p) of 
the initial product 2.2 is singular on the manifold n de- 
scribed by l g (p,0) = 0. The aggregate singular manifold of 
T K ° G(p) is U jeG K g . In general, the latter is singular in the 

sense of differential geometry, so one splits it into non-singular 
components Ky (labeled by an index y; each such component 

is a smooth open manifold). To each component there corre- 
sponds a subproduct Fy(p)cz G(p) (but unlike the Euclidean 

case, here different 7Cy may correspond to the same sub- 
product; e.g. the apex of the light cone and its two cones corre- 
spond to the same propagator). F Y (p) = y(p) contains all the 
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factors from G that are singular everywhere on Ky . Denote as 

G \y(p) the product of all factors from G that do not belong 
to y. One may say that y represents a "subgraph"; its corre- 
sponding product of propagators Fy(p) = y(p) is determined 

uniquely. 13 

Analytical structure of singularities 2.6 

Here we have to set rules for power counting and — simul- 
taneously — to define the so-called homogenization — a sec- 
ondary expansion needed in non-Euclidean situations to reduce 
the coefficients of counterterms of AO to power-and-log form. 
The discussion below is in the context of a given subgraph y. 
One considers a general point p Q of 7iy and the behavior of 

Fy ( p) when p — > p Q along directions that are transverse to 

7ty . The components of p that are tangential to 7ty are 

"spectators" and are to be ignored. After introducing appropri- 
ate local coordinates near p Q and redefinitions, we may as- 
sume in what follows that p does not have spectator compo- 
nents and that 7iy = {0} . 

Scaling Scale p i -> A Wi p i , ff — > A"* K so that: 

(i) all scaling exponents n a are positive integer; 

(ii) for each g £ 7 , l g (p,K) scales as 

A"s[/™ aln (/?,fc) + 0(A)], and the scaling exponent for K is the 
minimal value ensuring this for given n t ; 

(iii) if 7l™ ia is the manifold on which /™ in (p,0) = 0, then 
^main _ de f med as n ^main ; g me sgt of aU p such mat 

jmain^o) = for all g e y — coincides with n y = {0} . The 
latter means, in particular, that the collection of all 
l™ m (p,K) , g ey depends on all components of p. 

These properties ensure that each step of our expansion 
procedure is a mathematically correct transformation. A re- 
markable fact is that the scaling satisfying (i)-(iii) need not be 
unique (e.g. in the case of radiative corrections to the example 
of Sec. 3.1). The beautiful mathematical mechanism of how dif- 
ferent scalings result in the same final answer, is mentioned 
below after the definition of homogenization — we have al- 
ready had opportunities [15] to emphasize a remarkable stabil- 
ity of the method of AO that yields correct results even in 
counterintuitive cases as long as one applies it in a systematic 
manner. Different correct scalings do differ in the number of 
intermediate steps leading to the (same) final result. An opti- 
mal definition is as follows: Split I (p,0) — > L' (p) + L"(p) 

with L' g {p) = 0{p), L' g -(p) = 0(p 2 ). Then split 
p-^(X,Y,Z) so that all L' (p) depend on, and only on X 



and all l g (p,0) with L' g (p) = are independent of Z . Then X 

and Y scale with A 2 , and Z scales with A. 14 As a simple check, 
this rule correctly yields a uniform scaling in all components of 
p and k both when all l g (p,K) are linear functions, and when 

they are all purely quadratic functions (the Euclidean case). 

Power counting 15 One performs the power counting to de- 
termine the strength of singularity in each term of T^. °G(p) 

near generic points of each Uy . For that, one drops from de- 
nominators (and/or arguments of 5-functions) all but those 
components of l g (p,0) that scale with the lowest power of A, 

and introduces into the numerator the factor 
A 2dimX+2dimF+dimZ mat corresponds to the sca ii ng f the in- 
tegration measure. Collecting all powers of A one obtains an 
overall factor A my where C0y can be appropriately called the 
singularity index of the subgraph (e.g. (Qy =2 corresponds to 
quadratic divergence etc.). 

Homogenization 2.7 

The homogenization parameter E,y is introduced as fol- 
lows: (i) scale l g (p,K), gey as described; (ii) drop the over- 
all A" s ; (iii) replace A with <jL . The operation of homogeni- 
zation (denoted as H^, ) is as follows: (i) introduce ^ as just 
described; (ii) expand in E, y ; (iii) set S,y — » 1 in the result. 
Hy is meant to be applied to integrals similar to those we set 

out to expand from the very beginning (see below Eq.2.12). 
Therefore, the expansion in H^, must in general be performed 

in the sense of distributions — requiring the entire machinery 
of AO with another level of homogenization etc. Since at each 
step the expansion problem simplifies (dimensionality of the 
integration space is reduced), the recursion stops correctly after 
a finite number of levels of homogenization. In many interest- 
ing cases, however, non-trivial singularities requiring coun- 
terterms do not occur and this expansion degenerates into a 
simple Taylor expansion. 

The occurrence of secondary expansions with their corres- 
ponding homogenizations is behind the mechanism that en- 
sures independence of final results of the choice of scaling 
(provided the latter satisfies the three conditions mentioned in 
Sec. 2.6): non-trivial counterterms for secondary homogeniza- 
tions yield expressions that correspond to alternative scalings, 
whereas the "formal" part of the homogenization expansion 
yields zero by explicit integration. Demonstrating this mecha- 
nism in detail requires a much more detailed exposition than 
we can afford here. 



13 The pinch/non-pinch classification of singularities (the usual starting 
point of the conventional analyses; cf. e.g. [27]) is actually irrelevant for 
the analytical study of the singularities in general, and for the construction 
of AO in particular: the non-pinched singularities simply correspond to 
parts of singular manifolds where the corresponding counterterms nullify 
(which can be deduced e.g. directly from their expressions). Such a nullifi- 
cation is rather accidental from the point of view of analytical structure of 
singularities. 



14 For general non-quadratic (cubic etc.) functions one may recur to a 
technique based on Newton's polyhedra similar to that used in the theory 
of singularities of differential mappings [37]. 

15 This is regarded as a technical problem in the context of the theory of 
QCD factorization theorems [27], esp. in the case of non-leading power 
corrections and mixed soft/collinear singularities. Our rules should settle 
the issue (see also [30]). 
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Structure of asymptotic operation 2.8 

Recall the general formula for AO [15]: 
As K o G(p) = T„ o G(p) + £ r E r (p) x [T r ° G \ y(p)] . 2.9 

It is rather natural: each singularity of the formal expansion 
T^. °G(p) receives a counterterm E^(p) localized on the cor- 
responding singular manifold Ky . So, the summation runs 
over all subgraphs y, and Ey (p) have the form 

E Y (p) = Y, a E Y ,a(K)xSy,a(P), 2.10 

where summation runs over a complete set of 5-functions lo- 
calized on Ky , and Ey a (K) in general are — unlike the 

Euclidean case where each Ey a (K) is proportional to one in- 
teger power of K — series in (non-integer) powers of K with 
coefficients that are polynomials of logs: . As was already em- 
phasized, finding those coefficients is the central analytical 
task of the theory of asymptotic expansions of Feynman dia- 
grams. 

Consistency conditions for the coefficients 2.1 1 

The finding of Ey a (K) is performed according to the same 

general recipe as in [15], i.e. via consistency conditions. Sup- 
pose one wishes to construct the coefficients for E r a (K) for 

one subgraph T. One assumes that for all y < T the construc- 
tion has been performed (the usual ordering of subgraphs with 
respect to increasing codimensionality of singular manifolds is 
assumed here; cf. [15]). An appropriate choice of coordinates 
ensures that the singularity one is after is localized at p = . 
Then the coefficients are given by the following formulas (cf. 
sec. 12.3 in [15] and eq. (20.5) in [16]): 

E ra (K)= lim fdZdTdZ<5(ZA" 2 ,TA" 2 ,ZA" 1 ) 

x* a (p) H r [l-As^]or 

= Jdpff a (p)H { or(p,ic). 2.12 

(The polynomials v a form a complete dual set for the <S- 

functions 8 r a — exactly as in the Euclidean case. As'^. ° T 

differs from As^. ° T by absence of the term with y = T in the 

corresponding sum 2.9.) The first expression demonstrates how 
the intermediate cutoff is removed. Notice the asymmetry of 
the cutoff which corresponds to the asymmetry of the scaling. 
The role of He (a new element compared to the Euclidean 

case) is to split the coefficient into terms with pure power-and- 
log dependence on K (which motivates its definition). The sec- 
ond expression takes into account the nullification of the sub- 
tracted ("shadow") terms in dimensional/analytic regulariza- 
tion, which is due to their pure-power behavior under the 
(asymmetric) scaling. 16 ' 17 The power of K for each term in the 



16 Nullification of integrals of the form Jd°p p 2a = is well-known. It 

is due to the fact that the dimensional regularization preserves formal 
scalings. In the new situation, we are dealing with zero integrals of, very 

roughly, the form J d D ~ 2 p±dp + + pY) =0. The underlying 
reason for their nullification is the same — preservation of scaling proper- 



last expression in 2.12 is determined by scaling out K using the 
scaling rules for r already explained. 

An important point to remember is that the construction of 
asymptotic expansions (including AO) is always, strictly 
speaking, carried out for a particular finite precision 0(K N ) . 
This implies that one subtracts only the series to that precision 
in the integrand of 2.12. Correspondingly, the counterterm is 
also defined at that step only within that precision, whereas 
expansion implied by He — if carried too far — would gener- 
ate terms of excessive precision 0(K N+n ) — terms that would 
also be divergent in the UV region! One can verify by power 
counting, however, that all the contributions of precision 
0{k n ) are exactly those whose (formal) UV convergence is 
ensured by the subtraction. In the end one can forget about the 
restriction 0(k n ) because the formula (the last expression in 
2.12) is independent of it. 18 

Another important point concerns the diagrammatic inter- 
pretation of the subtraction in the first expression in 2.12. In 
the Euclidean case, it was shown [34; 21] that the subtractions 
of this sort exactly correspond to the standard Bogoliubov UV 
i?-operation — in agreement with the fact that the integrals 
that occurred there for E T a were exactly diagrams with local 

operator insertions. In the general non-Euclidean situation, 
there is no universal (i.e. valid for all regimes) operator char- 
acterization for the integrals that occur after the homogeniza- 
tion (the expansion He distorts standard propagators in differ- 
ent ways for different asymptotic regimes). But the prescription 
for their UV renormalization is always determined uniquely by 
the structure of subtractions in 2.12. 

As a last remark, an interesting technical point may be 
mentioned. Namely, despite the presence of the additional ex- 
pansion He in Eq. 2.12 compared to the Euclidean case, when 

one integrates out the 5-functions similarly to the procedure of 
sec. 5.3 of [8] to obtain the AO in integral form, this additional 
series blends into a single Taylor expansion with the expansion 
resulting from derivatives of 5-functions (cf. eq.(5.8) in [8]). 
This is easily explained if one notices that the non-Euclidean 
prescriptions remain valid in the general non-linear/non- 
uniform case, so that the "spectator" product G\T can be 
formally absorbed into T and treated as a "deformation" of the 
latter, and when obtaining the AO in integral form, only those 
counterterms of the AO would survive that are proportional to 
5-functions without derivatives. But then the expansion He 

would be the only expansion remaining in the end! (Note, 
however, that the expansion implied by He affects the entire 

G — in contrast with the Euclidean case where only certain 
subgraphs are thus expanded; cf. the expansions of "heavy 



ties by the dimensional/analytical regularization, although now the scaling 
is non-uniform in different components of p . 

17 To obtain the corresponding regularization-independent formula one 
replaces ASk — > ty- ° ASk where tf is an appropriate generalization of 
the corresponding Euclidean subtraction operator; see [16] and [30]. 

18 However, a great caution must be exercised when establishing corre- 
spondence between the terms of the last expression in 2.12 and the singu- 
larities of the formal expansion T K »G. The correspondence is rather 
tricky which fact caused some confusion in the literature (see a discussion 
in [19]). 
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knots" in [8].) An amusing implication is that use of the gen- 
eral non-Euclidean results would actually simplify the con- 
struction of Euclidean AO in integral form. 

Examples 3 

The two examples we consider are related to well-known 
simple integrals so that checks are possible, yet involve ex- 
plicit phase space so that a systematic treatment via any other 
method would be problematic. In both examples, obtaining 
higher terms of the expansion is tedious but entirely straight- 
forward — an exercise that is left to an interested reader. Sub- 
stantially more involved examples will be presented elsewhere 
[38; 39]. 

Threshold Q 2 ~Am 2 for the kinematics Q^m+m 3.1 

This example corresponds to a near-threshold creation of a 
pair (e.g. e + e~ by a photon). Remember that asymptotic op- 
eration commutes with multiplication by polynomials so that 
non-scalar numerators are ignored in our example. 19 The phase 
space is represented as 



in Taylor-expanding with respect to the quadratic terms in p<j 
etc.), we find: 



\d D p w(p) S + (p 2 - m 2 )S + ((Q - p) 2 - m 2 ) . 



3.2 



where 5 + (k —m ) = 0(k )S(k —m ) , Q is the momentum 
of the "photon" that decays into the pair, and the arbitrary 
weight w (corresponding to arbitrary cuts experimentalists may 
use for event selection) means that the phase space is treated 
as totally exclusive, so that one essentially deals with the prod- 
uct of O's and <5's interpreted as a distribution. We are going to 
extract the first non-trivial (square root) contribution using the 
prescriptions set forth above. 

To describe the asymptotic regime we choose 
Q = (2m + K,0) with fc^O.Then 

K = (Q 2 - Am 2 )/4m - {Q 2 - Am 2 ) 2 / (Am) 3 + ... 

It is convenient to perform a shift p — > p + in, rh = (m,0) . The 
formal expansion of the product in K then is: 

<S + ((p + m) 2 -m 2 ^5 + ((m + K - p) 2 -m 2 j 

= 8 + {pl + 2mp Q - p 2 )8 + (p% - 2mp Q - p 2 ) 

+2>r[-p + m] 8 + {p 2 + 2mp Q - p 2 ) 8' + {p 2 - 2mp Q - p 2 ) +. . . 3.4 

The singularity at the point p = is seen from the fact that the 
arguments of the two 5-functions reduce to the form 
S(p Q )S(p 2 ) , with the second <5 ill-defined in the case of 

higher derivatives (the singularity is regulated by the dimen- 
sional regularization). To do the power counting according to 
the prescriptions of Sec. 2, one identifies X <-> (p ,K\), Y <-> p 

with Z empty. The first product of 5's is convergent whereas 
the second one is linearly divergent. In any case, a counterterm 
of the form E(K)8(p) is necessary. Evaluating E(k) ac- 
cording to the given recipe (the homogenization consists simply 



Of course, they do play a role in determining which particular terms 
from the asymptotic expansion contribute to a particular process at a par- 
ticular precision with respect to the small parameter. But our aim here is 
only to demonstrate the essential analytical mechanism. 



E(K) = \& D p8{2mp - p 2 )s{2m(K - p Q ) - p 2 ) 



_ 1 v K 



1/2-6 



- . r x(m/c) 

4m r(3/2-e) 

where e = 4(4- D) . The final result for 3.2 is as follows: 



3.5 



jd D pw(p)S + (p 2 -m 2 )8 + ((Q-p) 2 -m 2 ) 

= -f^^Q 2 ~ 4m2 w(in) + o{Q 2 -Am 2 ). 3.6 

(The integral of the first term of the formal expansion — the 
first term on the r.h.s. of 3.4 — is zero for smooth w.) This re- 
sult is checked by noticing that the case w = 1 corresponds to 
the imaginary part of a simple self-energy diagram; cf. the ex- 
plicit result in [40], eq.(24.5). 



The behavior of 2->2 at s->-h=o 



3.7 



Our second example corresponds to the matrix element 
squared with exclusive phase space for k + + k_ — > k' + + k'_ 

(where k' ± = k ± ± p ) via simplest /-channel exchange of a 
scalar particle with momentum p at large s = (k + + k_) 2 . The 

example is motivated by the large- y/small-x problem in QCD 
[41], so all particles (partons) are massless and the external 
partons are slightly off-shell (time-like). The expression to be 
expanded is: 

s\&p w(p)8 + [(k + + p) 2 )8 + [(k_ - p) 2 )p- 2 , 3.8 

where <5 + ((fc + + p) 2 ) = 0(k° + p°)S((k + + p) 2 ) etc., and we 

have replaced for simplicity the exchange particle's squared 
propagator p~ 4 by p~ 2 (imagine e.g. that there occurred a 
3.3 cancellation with p 2 from a non-scalar numerator), w is an ar- 
bitrary weight that describes cuts, observables etc. The as- 
ymptotic regime is described by s —> °° with k+, p 2 = 0(1) . 
Introduce two light-like vectors q ± such that 2q + q_ = 1 and 

k ± = *Js q+ + 0(s~ 1 ' 2 ) . Taking out *J~s from each 5 + and de- 
noting k = s~ 112 , rewrite 3.8 as 
j dp w( p) p~ 2 8 + (2q + p +Kkl + Kp 2 + 0(K 2 p) j 

8 + (-2q_p + i& 2 +Kp 2 + 0(K 2 p)Y 3.9 
The first term of the formal expansion is: 

\dpw(p)8 + (2q + p)8 + (-2q_p)p- 2 + 0(K) . 3.10 

Singularities are localized on n ± = {p = zq±\ z e D&} and 

7t = {0} . The two collinear singularities Jt ± are similar to 

those considered in [23]. Following the above rules, one finds 
the appropriate scaling (in the Sudakov coordinates): 

Pp — > A 2 p+ , pj_ — > Ap x , k — > ??K (for n + the component 

Pp is spectator and does not scale). One finds that in the 

leading order the two singularities are logarithmic, so only 5- 
functions without derivatives have to be added to 3.10: 

dzcAz,K)8(p-zq + )+\ d Z c_(z,K)8(p- zq_) . 3.11 



F. V. Tkachov Near-threshold expansions Phys. Lett. B412 (1997) 350 [hep-ph/9703424] 



7 



The corresponding coefficients are easily found: 
C ± (Z,K) = jd D p8 + (±2pq ± +Kkl)8(z-2pq+)p- 2 

~ iW*Kf ■ 3.12 

where A" = G(A > 0)A a and we have dropped irrelevant coef- 
ficients. However, when one multiplies the collinear coun- 
terterms by the corresponding G \ y+ = 5 + (+2q T p) , one finds 

that the product vanishes because c ± (z) ~ z~ e (cf. the case 
with a non-zero mass where a similar product was ill-defined 
[23]). 

The scaling for the singularity localized at Jt Q = {0} is uni- 
form in all components p — > hp , K — > Xk . The counterterm 
then is c Q (K)S(p) with 

C (K) = Jd D p 5 + {2q + p + Kkiy + (-2q_p + Kfc 2 ) p~ 2 

~jK~ 2£ . 3.13 

Due to all the <5-functions, the integrals in this example are 
performed very easily. 

Finally, the expansion is: 

Eq. 3.8 = Eq. 3.10 +c (k)w(0). 3.14 

Performing integration in 3.10 with w=l one finds that the 
poles e -1 in 3.14 cancel (as expected). One sees that the 
log/c contribution is associated with the diffractive region 
p~0. 

If one had normal non-cut propagators instead of the two <S- 
functions in 3.8, the vanishing of collinear counterterms would 
not occur, which would give rise to the well-known log 2 k 
terms (taking imaginary part eliminates one logarithm). Oth- 
erwise the integrals are just a bit more involved, and the cor- 
rectness of calculations is easily checked by comparison with 
the well-known analytical results (see e.g. [42]). 



For Monte-Carlo type applications one needs to recast the 
results into a regularization-independent form. The correspond- 
ing prescriptions are derivative from the power-counting rules 
of Sec. 2.6 and follow the general pattern of [16; 21]. More 
details will be given elsewhere [30]. 

Conclusions 4 

The recipe for the non-Euclidean (near-threshold) asymp- 
totic operation given above is a very general one, is valid prac- 
tically for any non-Euclidean asymptotic regime, and for any 
individual Feynman diagram, including unitarity diagrams with 
cut propagators as well as diagrams in non-covariant gauges, 
heavy-quark and non-relativistic effective theories, etc. But de- 
spite its generality, it relies on few analytical principles, which 
is important for calculationists who wish to have a complete 
intellectual control over what they are doing (cf. e.g. 
[22; 24; 26]). In many applications (namely, automated higher- 
order calculations and derivations of all-order near-threshold 
expansions in operator/effective Lagrangian form) one would 
need the rules for AO in integral form similar to the Euclidean 
ones found in [5-8], that are currently in wide use (cf. [9-11]). 
Such rules are not difficult to obtain — but then the universal- 
ity of formulas is lost and each regime has to be treated sepa- 
rately. Studying the operator form of expansions for various 
non-Euclidean (near-threshold) regimes and performing the 
corresponding calculations will be a topic of active research in 
the coming years [26; 38; 39; 43]. 

Acknowledgements I am indebted to V.I. Borodulin, G. Jikia 
and N. A. Sveshnikov for an extremely important support. 
I thank D. Broadhurst and P. Nason for pointing out the impor- 
tance of studying non-zero thresholds in QED and the Standard 
Model, S.Jadach for emphasizing the importance of exclusive 
treatment of phase space, A.V. Radyushkin for bibliographic 
advice, and F. Berends for an encouraging laugh. This work 
was supported in part by the Russian Foundation for Basic 
Research under grant 95-02-05794. 



F. V. Tkachov Near-threshold expansions Phys. Lett. B412 (1997) 350 [hep-ph/9703424] 



8 



References 

[1] F.V. Tkachov, in Quarks-82. Proc. Int. Seminar. Sukhumi, USSR. 

May 5-7, 1982. INR, USSR Acad. Sci. 
[2] F.V. Tkachov, Sov. J. Nucl. Phys. 56 (1993) 180 [hep- 

ph/9701272]. 
[3] F.V. Tkachov, Phys. Lett. 124B (1983) 212. 
[4] S.G. Gorishny, S.A. Larin, and F.V. Tkachov, Phys. Lett. 124B 

(1983) 217. 

[5] G.B. Pivovarov and F.V. Tkachov, preprint INR P-370 (1984) 
[KEK no. 8502210]. 

[6] G.B. Pivovarov and F.V. Tkachov, preprint INR Fl-459 (1986) 
[KEK no. 8610360]. 

[7] G.B. Pivovarov and F.V. Tkachov, in Quarks-86. Proc. Int. Semi- 
nar. Tbilisi, USSR. May 1986. VNU Science Press, 1987. 

[8] G.B. Pivovarov and F.V. Tkachov, Int. J. Mod. Phys. A8 (1993) 
2241 [hep-ph/96 12287]. 

[9] S.A. Larin, F.V. Tkachov, and J.A.M. Vermaseren, Phys. Rev. 
Lett. 66(1991) 862. 
[10] S.A. Larin, F.V. Tkachov, and J.A.M. Vermaseren, Phys. Lett. 

272B (1991) 121. 
[11] J.H. Kiihn, in Cracow International Symposium on Radiative Cor- 
rections, Cracow, Poland, 1996 [hep-ph/96 12372]. 
[12] F.V. Tkachov, Int. J. Mod. Phys. A (1998, in print) [hep- 
ph/9703423]. 

[13] N.N. Bogoliubov, Doklady USSR Acad. Sci. 82 (1952) 217. 

[14] N.N. Bogoliubov and DV. Shirkov, Introduction to the Theory of 
Quantized Fields. Interscience, 1959. 

[15] F.V. Tkachov, Int. J. Mod. Phys. A8 (1993) 2047 [hep- 
ph/96 12284]. 

[16] A.N. Kuznetsov, F.V. Tkachov, and V.V. Vlasov, preprint PSU- 

92/108-T (1992) [hep-th/96 12037]. 
[17] R.D. Richtmyer, Principles of Advanced Mathematical Physics, 

vol. 1. Springer- Verlag, 1978. 
[18] F.V. Tkachov, preprint INR P-332 (1984). 

[19] F.V. Tkachov, in The VIII Int. Workshop on High Energy Physics. 
Zvenigorod, Russia. 15-21 September 1993. INPMSU, 1994 [hep- 
ph/96 12288]. 

[20] F.V. Tkachov, Phys. Lett. 125B (1983) 85. 

[21] A.N. Kuznetsov and F.V. Tkachov, preprint INR-809/93 (1992) 
[hep-th/9612038]. 



[22] V.P. Spiridonov and F.V. Tkachov, Phys. Lett. 260B (1991) 109. 
[23] J.C. Collins and F.V. Tkachov, Phys. Lett. 294B (1992) 403. 
[24] V.I. Borodulin and G. Jikia, preprint Freiburg-THEP 96/19 

(1996). [hep-ph/9609447] 
[25] J.C. Collins, L. Frankfurt, and M. Strikman, preprint CERN- 

TH/96-314, PSU/TH/168 (1996) [hep-ph/9611433, version of 20 

Dec 96]. 

[26] V.I. Borodulin and G. Jikia (to be published). 

[27] G. Sterman, preprint ITP-SB-96-4 (1996) [hep-ph/9606312]. 

[28] NA. Sveshnikov and F.V. Tkachov, Phys. Lett. B382 (1996) 403 

[hep-ph/95 12370]. 
[29] F.V. Tkachov (in preparation). 
[30] F.V. Tkachov (in preparation). 

[31] G. 't Hooft and M. Veltman, Nucl. Phys. B44 (1972) 189. 
[32] E.R. Speer and M.J. Westwater, Ann. Inst. H.Poincare A-XIV 
(1971) 1. 

[33] A.N. Kuznetsov and F.V. Tkachov, talk at Quarks-88. Proc. Int. 

Seminar. Tbilisi, USSR. May 1988. 
[34] A.N. Kuznetsov and F.V. Tkachov, in Renormalization Group '9 1 . 

II Int. Conf. Dubna, USSR. 3-6 Sept. 1991. World Scientific, 

1992. 

[35] G. 't Hooft, Nucl. Phys. B61 (1973) 455. 
[36] F.V. Tkachov (in preparation). 

[37] V.I. Arnold, A.N. Varchenko, and S.M. Guseyn-Zade, Singularities 
of differentiable mappings. Classification of critical points, caustics 
and wave fronts, vol. 1. NAUKA, 1982 (in Russian). 

[38] P.S. Cherzor, NA. Sveshnikov, and F.V. Tkachov (in preparation). 

[39] N.A. Sveshnikov, F.V. Tkachov, and N.I. Ussyukina, (in prepara- 
tion). 

[40] N.N. Bogoliubov and D.V. Shirkov, Quantum Fields. NAUKA, 

1980 (in Russian). 
[41] L.N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 642. 
[42] N.I. Ussyukina and A.I. Davydychev, Phys. Lett. 298B ( 1 992) 

363. 

[43] A.V. Stepanova and F.V. Tkachov , 'Algebraic algorithm for next- 
to-leading order calculations in the large-s/small-x regime", in: 
Proc. Int. Workshop QFTHEP'97, Samara, Sept. 4-10, 1997 [hep- 
ph/97 10242]. 



